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Abstract. In this article we provide numerical and analytical evidence that 
some degenerate dispersive partial differential equations are ill-posed. Specifi- 
cally we study the K(2, 2) equation ut = (u 2 ) xxx + (u 2 ) x and the "degenerate 
Airy" equation ut = 2uu xxx . For K(2, 2) our results are computational in na- 
ture: we conduct a series of numerical simulations which demonstrate that data 
which is very small in H 2 can be of unit size at a fixed time which is indepen- 
dent of the data's size. For the degenerate Airy equation, our results are fully 
rigorous: we prove the existence of a compactly supported self-similar solution 
which, when combined with certain scaling invariances, implies ill-posedness 
(also in H 2 ). 



1. Introduction 

Dispersion plays a pivotal role in the analysis of the governing equations for a 
large number of physical phenomena, ranging from the evolution of surface water 
waves to the formation of Bosc-Einstcin condensates. A common theme in such 
problems is that uniform dispersive effects frequently control nonlinear terms which 
may otherwise lead to ill-posedness. Consider for instance the following class of 
quasilinear equations of KdV type studied in [8J : 

(1) u t = 93U XXX + 92U xx giu x ~t~ gou -\- h 

where gj = gj(d%.u, u) and h = h(x, t). There it is shown that such equations are 
well-posed provided |g 3 | > c > and g 2 > 0. The first of these conditions guaran- 
tees that dispersive effects from g-iU xxx do not vanish while the second prohibits the 
term g 2 u xx from acting like a backwards heat operator. Similarly, the results on 
existence of solutions to the quasilinear dispersive equations found in [T5] (about 
Schrodinger equations), [T] (water waves) and [25] (magma dynamics) share the 
feature that the dispersive effects arc nondegenerate, to name just a few examples. 

On the other hand, there are a number of physical problems in which the mecha- 
nism which causes the dispersive effects is not only nonlinear but degenerate. Some 
examples can be found in the study of sedimentation [25], shallow water 0, gran- 
ular media [TH] HH [2U] , numerical analysis [2] and elastic rods and sheets H5J . 
While there are a large number of articles in which degenerate dispersive equations 
(henceforth referred to as DDE) are derived [H2 [2Ql EJ Mi U ES] , special solu- 
tions computed [23] [2SJ [TTJ [25] [T5] or numerical computation of solutions performed 
[251 [271 [211 [TU1 [T71 [5J , the existence theory for the initial value problem of these 
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equations is largely undeveloped^] The articles [TOl 02 [13] prove some a priori es- 
timates for DDE. [2] and [2] prove existence (but not uniqueness) of solutions to a 
very special class of DDE. [T2] contains a semi-rigorous discussion of existence issues 
for these sorts of equations, including defining "5-entropy solutions" by borrowing 
ideas from first order conservation laws. None of these articles fully settles the 
well-posedness issues for degenerate dispersion, and thus our interest is in studying 
the Cauchy problem for some simple and common DDE. 
In particular we study the following two equations: 

(2) u t = 2uu xxx + 6u x u xx + 2uu x 
and 

(3) u t = 2uu xxx 

for i£R. The first of these is the K{2, 2) equation of Rosenau & Hyman [24] . Their 
goal was to understand the role of degenerate nonlinear dispersion in the formation 
of patterns and this equation (together with other members of the K{m, n) family 
Ut — {u m ) x + (u n ) xxx ) has come to serve as the paradigm for degenerate dispersive 
systems. The most celebrated feature of this equation is the existence of compactly 
supported traveling waves, a.k.a. "compactons." These are given by 

, n ^ / , n 4A 9 fx + Xt\ , , . 
u(x, t) = Q x (x + Xt) := — cos^ I - j X i(x + Xt) 

where xi 1S the indicator function for the interval / = [— 2ir, 2ir] and A ^ 0. More 
recently, K(2,2) has been used to model pulses in both elastic and granular me- 
dia [El [21] . Equation ^ , a "degenerate Airy equation" , appeared in [13] and is 
amongst the simplest equations which could be said to feature degenerate disper- 
sion. We study it here due to its similarity to the K(2, 2) equation and because it 
is of the form of ([!]) but docs not meet the uniform dispersion hypothesis needed 

Surprisingly, our conclusions are that both equations are ill-posed for data in H 2 . 
In particular, their solutions do not depend smoothly on their initial conditions. 
Our results for J2j) are computational in nature: we conduct a series of numerical 
simulations of (12^ which demonstrate that data which is very small in H 2 can be 
of unit size at a fixed time which is independent of the data's size. For (|| our 
results are fully rigorous: we prove the existence of a compactly supported self- 
similar solution which, when combined with certain scaling invariances, implies ill— 
posedness. Note that our focus on H 2 initial data is due to the fact the compacton 
solutions of -&T(2,2) have a discontinuity in their second derivative at the edge of 
their support [TH]- Consequently, they reside in H 2 but not H 3 , and so it is natural 
to work in this space. 

Our suspicion that ([2| is ill-posed has its origin in the various numerical sim- 
ulations of collisions between compactons in K(m,n) equations. A wide variety 
of numerical methods has been used to compute solutions: pseudo-spectral [2~4"] . 
particle [B], Pade [55J [37], discontinuous Galerkin [T7], and finite difference [TU] , 
Each of these simulations shows that the computed collisions between compactons 



^The notable exception to this is the well-studied Camassa-Holm equation ut — u xx t = —3uu x + 
2u x Uxx + uu xxx which when rewritten as a nonlocal evolution equation could be said to be 
degenerate and dispersive. This equation is integrable and consequently many results for it do 
not obviously generalize. 
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are nearly elastic: the compactons emerge from the collision nearly identical to 
their original state but phase shifted. As is common in collisions between soli- 
tary waves in "nearly" integrable PDE, the simulations show that in addition to 
the outgoing waves there is a small amplitude disturbance behind the waves. In 
nondegenerate problems, similar disturbances are sometimes called "dispersive rip- 
ples" and are viewed as being linear effects — the tail (approximately) solves the 
linearized equation [301 [3T] . However, since the linearization of |2| about zero is 
trivial this heuristic cannot apply. The disturbance is described variously as: a 
compacton/anti-compacton pair in [24] ; a numerical artifact in [281 110] ; and as a 
shock in [27]. No such disturbance is present at all in |B]. Lastly, the authors of 
[17] describe a "high-frequency oscillation" but go on to state "we do not know 
what is the source of these oscillations..." In each case, the qualitative features of 
the disturbance are quite different. Our feeling is that these discrepancies are not 
related to differences in the numerical methods or the implementation thereof but 
rather are due to an underlying issue with the equation itself. 

Moreover, inspection of the right hand side of ^ provides evidence that the 
equation may be ill-posed. Observe that when u x < 0, the term 6u x u xx is formally 
a backwards heat operator. The key question, initially raised in [24] . is this: can 
dispersive effects due to 2uu xxx counteract the instability caused by the backwards 
heat operator 6u x u xx ? Thus we should seek to understand the dispersive effects 
due to 2uu xxx — hence our interest in ([3]). Our results show that the answer to the 
key question is "no" . As we stated above, our results show that 2uu xxx is, in its 
own right, just as problematic as the backwards heat term. Differentiation of ^ 
with respect to x yields for w := u x 

w t — 2uw xxx + 2u x w xx . 

Note that the second term on the RHS above is also a backwards heat operator 
when u x < 0. The backwards heat term evident in ([2| is, in this way, hidden 
in ([3]). Thus we expect blow up of u x when u x is negative. Moreover, this blow 
up should be particularly catastrophic when the function crosses the x-axis with 
negative slope. 

This motivates our approach: the self-similar solution constructed has such a 
transverse crossing and the resulting ill-posedness is due to blow up in the first 
and second derivatives of the solution. While our proof that ^ is ill-posed does 
not naturally carry over to ([2]), the numerics are highly suggestive. Most tellingly, 
since the numerics are performed using strictly positive initial data, they suggest 
the ill-posedness can manifest even in the absence of such a transverse crossing. 

Remark 1. In 7\ 7 Craig & Goodman demonstrate that the differential equation 
u t = —xu xxx is ill-posed, whereas u t = xu xxx is not only well-posed but the solution 
map is essentially infinitely smoothing. The proofs of these facts follow from explicit 
formula for solutions of the initial value problems which can be constructed for each 
equation by mean of the Fourier transform and the method of characteristics. 

The remainder of this paper is organized as follows. In Section[2]we show that the 
existence of a self-similar solution to ^ gives, as a consequence, the ill-posedness 
of the equation. In Section [3] we construct the self-similar solution using tech- 
niques from spatial dynamics. In Section [4] we discuss our numerical strategy for 
demonstrating ill-posedness of p]) and present the results of our simulations. 
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2. Self-similar solutions imply ill-posedness for ^ 
We will prove there is a self-similar solution for ([3| of the form 

u(M) = A ((i4)V3, 

where A gives the profile. Though other choices for the scaling are possible, we 
choose this particular set so that our proposed solution exhibits blow up (as t — > 1~) 
in the first and second derivative while the L°° norm remains bounded. Inserting 
this into ^ yields: 

(4) 

Here r = x(l — i) -1 / 3 is the self-similar coordinate upon which A depends and we 
have (□)' = d/dr(n). We specify A(0) = and A'(0) = -1 so that the self-similar 
solution crosses the x-axis transversely with negative slope, the "bad" situation 
identified in the Introduction. 

We prove the following proposition concerning solutions of Q: 

Proposition 1. There exists /x* > 0, > 0, C* > 1 and a function A*(t) G 
C 3 ([0, r*)) with the following properties: 

(i) A*(t) solves Q for all re [0,r*). 

(ii) A*(0) = 0, A* (Q) = -1 and A'l(0) = /i*. 

(hi) A„(t) is negative in (0,t*) and has a unique minimum in the interior of this 
interval. 

(iv) lim_ A* (t) = lim_ A'„ (r) = 0. 

r— fr* t — vt^ 

(v) C^ 1 | log \t - r*|| < |A'„'(r)| < C* |log |r - t*\\ for all r € (0,r») sufficiently 
close to . 

The proof of this proposition is technical and difficult and we postpone it until 
the next section. First we utilize it to prove our main result for ([3|: 

Theorem 1. The Cauchy problem for ([3]), posed in H 2 , does not depend smoothly 
on the initial data. 

Proof. (Theorem [I]) The proof of the theorem is classical: we construct a sequence 
of solutions which are arbitrarily small in H 2 at t = but which are arbitrarily 
large when t = 1. Let 

(A*(t) ifre[0,n) 
A(t) := I —A*(—t) if r e (-t* , 0] and 
[ otherwise. 
This is the profile function for our self-similar solution. See Figure [T] 
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FIGURE 1. Numerically computed profile A(t). Note \i = ^ ~ 0.354875. 



Note that it is compactly supported, not unlike the famous Barenblatt solution 
of the porous medium equation [3] . Also observe that A" is discontinuous at r = 
and diverges logarithmically at ±t*. Nevertheless, these singularities are sufficiently 
mild so that A £ H 2 . Also observe that tA' is continuous for all r and consequently 
so is AA'" . Thus the function 

u±,i(x, t) := A 



(l-i)V3 

is in fact a classical solution of pi). That is to say, while d 3 u\^\ ^ L 2 , ui^d 3 u\.\ is. 
Compare this with the compacton solution of K(2,2): it is not in H 3 , though its 
square is. 

For any T € R and A > 0, (f3l) is invariant under the following transformation 
(5) U ^ t )^\u(j^,t-T). 

Thus we have a two-parameter family of self-similar solutions of the form: 

x 

AV3(T - 1) 1 / 3 



u\,t(x, t) := XA 
An elementary computation shows that 
(6) |KtM)II^ =A 7 /6( T -t) 1 /6|| A || L2 

+ A 5 / 6 (T - t)-^ 6 \\A'\\ L 2 + A J / 2 (T - t)- 1 / 2 !!/!"!!^. 
Fixing < e « 1, setting A(e) = e 2 and T(e) = 1 + e 16 results in 

C 

ll"A(e),T( e )(- ; 0)\\ H 2 < Ce and ||uA(e),T( £ )(-, 1)||h 2 > — • 
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Figure 2. Numerically computed solutions of ([7]) for increasing 
values of fi. 

(Here C > is a nonessential constant.) This pair of inequalities demonstrates that 
we have constructed a sequence of initial data for ([3| which is arbitrarily small but 
whose solution at time t = 1 is arbitrarily large. Thus the equation is ill-posed in 
H 2 . 

□ 

Remark 2. In |13j . the authors investigate self-similar solutions of the DDE 
u t = (uu x ) xx . In particular they numerically compute profiles for shock and rar- 
efaction solutions which are, in many respects, quite similar to the ones we find 
here. Namely, the shock solutions cross the x-axis with negative slope and form a 
singularity in u x at the time of the shock formation. The solutions they find are 
not in L 2 (R) — either they do not converge to zero as x — > oo or their decay is not 
rapid enough to be integrable. Our focus here is on well-posedness issues in H 2 and 
consequently our strategies and methods will differ from those of |13] . 

3. Construction of the self-similar solution 

This section is devoted to proving Proposition [I] We rewrite Q as the initial 
value problem (I VP): 

(7) &AA'" = tA' with A(0) = 0, A'(0) = -1, A"(0) = /x. 

Our strategy will be as follows: first we will show that the IVP Q has a unique 
solution, which depends on [i continuously. Next we will show that for sufficiently 
large the solution of Q goes to zero in finite time. Then we will show that for fx = 
there is a smallest time < r 2 < oo such that the solution of Q has its first local 
maximum at r 2 and remains strictly negative for all r e (0, r 2 ]. See Figure[2j Then 
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an open/closed argument will give the special solution A* for a choice fi = /i* > 0. 

The differential equation Q is nonlinear, nonautonomous and singular at A = 
0. As such, carrying out the details of our strategy is quite challenging. The 
key observation for proving Proposition [l] is that Q is, in a loose sense, close 
to being hamiltonian. Specifically, one can rewrite the right hand side of p} as 
t d 

In Ij4(t)|. If we then consider the related equation 

6 dr 

(8) B>" = ^\n\B(r)\ 

obtained from Q by replacing the nonautonomous r/6 with the constant 1, one 
sees that this equation is in fact hamiltonian. Specifically one can integrate both 
sides to find: 

s " (r) = Z^ ((B " (0) ~ 1)B + i31n|i?l) 

The conserved energy is 

~(B'f-((B"(0)-l)B + Bln\B\) 

and most of the dynamics of ^ can be determined by analyzing this quantity. 
In the following proof of Proposition [l] we will frequently make use of quantities 
similar to this energy which, though not conserved for solutions of Q, do provide 
quite a bit of dynamical information about solutions. 

3.1. Well-posedness of |7]). Now we prove ^ has unique solutions and that such 
solutions depend smoothly upon \i. 

Lemma 2. For any /i 6 R, there exist r > and S > such that for each 
/i G (fJ,Q — S, /io + 8) the IVP |7| has a unique solution A : [0, To] — > R, which 
satisfies A(t) < for all r G (0, To] . In addition, there exists C > such that for 
any fii,fi2 G (/io — 5, Ho + 8) the solutions A\{t) and ^(t) to the IVP Q with 
A'{(0) = [i\ and A'^O) — [i2> respectively, satisfy 

\A 2 (r) - A!(t)| + \A' 2 (r) - A[{r)\ + \A»(t) - A'{(r)\ < C\^ 2 Ml | 

for all t G [0, To] • 

tA' 

Proof. Letting F(A, A', A", t) := — , Q can be written as A'" = F(A, A', A", t). 

If F(A, A' , A" ,t) were continuous in t and Lipschitz in (A, A', A") throughout 
a neighborhood of (A, A', A", t) = (0,-1,^,0), then Lemma [2] would be a triv- 
ial consequence of the standard theory on existence and uniqueness of solutions 
and their continuous dependence on initial conditions. Unfortunately, in our case 
F(A, A' , A" , t) is not even defined if A = 0. However, notice that if there exists a 
solution A(t) to the IVP Q then L'Hopital's rule implies: 

1-^ = -- 
r->o+ 6A(t) 6 

Therefore our strategy will be to remove the singularity at (A, t) = (0, 0) by modify- 
ing F in such a way that the resulting modification is continuous in a neighborhood 
of the initial data and agrees with F for solutions of §f§. 
Fix /io G R and k > 0, and let 

p(r) := -t + ^(ju - k)r 2 and q(r) := -t + ^(/i + k)T 2 . 
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Given the initial conditions in |7]), p(r) and q(r) formally give crude lower and 
upper bounds on A(r) for r close to 0. Now construct the following function: 



tA' 

6p(r) 
tA> 

F{A,A',A",t) :=\% 

6qV) 
_A!_ 

IT 



if t > and A < p(r), 

if t > and p(r) < A < q(r) 
if t > and q(r) < A, 

if t < 0. 



F is well-defined and continuous throughout a small neighborhood of (A, A' , A" , r) = 
(0, — 1, /Lt, 0) for any /i6R, Then we have the existence of a solution A(t) to the 
I VP: 

A'" = F(A, A', A", t) with A(0) = 0, A'(0) = -1, A"(0) = ft. 

Furthermore, it is easy to verify that there exist t > and < 6 < k such that for 
each /i G (/i — ^ A*o + <5) the solution A(t) is defined for all r 6 [0, r ] and satisfies 
p(/r) < A(t) < q(r) < for all t G (0, r ]. It follows that the restriction of A{t) to 
[0, To] is in fact a solution to the IVP (FtI). 

Note that along any solution of Q , ^-^ ( J -^ooasr^0 + . Consequently, 

it is not possible to modify Ftoa function that is Lipschitz in A in a neighborhood 
of (A, A' , A" ,t) = (0, — 1,(4, 0). Therefore, we have to establish the uniqueness of 
the solution by other means, namely the nearly hamiltonian structure of Q. 

For any fii,fi2 G (Mo — fio + S) , let A\ (r) and A 2 (r) be solutions to the IVP (J7]) 
with A/(0) = /Ux and A 2 '(0) = /i 2 , respectively. Recall that our choice of 8 and r 
ensures p(r) < Ai, 2 (t) < q(r) < for all r G (0, To]. Thus we can take a sufficiently 

small T G (0, min{l, To}], independent of pi and /i 2 , such that j4i. 2 (t) < — — r < 
for all r G (0, T] and 



Then for any t G (0,T], 
(9) 



M 



In 



max 

t£(0,T] 



A 2 (r) 



Mr) 



Mr) 



< 



AI 



1 - M 



< 



< 2M, 



.Mr) 

where the equalities hold if M = 0. 

Since both Ai and A 2 are solutions to the ODE, we have 



A!£ - A'" = 



A', 



6 \A 2 

Integrating the above equation gives 



-A 

Ax 



Ax 
A 2 



M 
Ai 



Mr) -Mr) 



6 



In 



A 2 (r) 
Ai(r) 



In 



Ms) 
Ms) 



ds + 



where S/i := /i 2 — fi\. Then by incorporating ([9]), we obtain that for all r G (0,T], 



(10) 



^(r)-A'/(r)|<-Mr 
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Together with the initial conditions Ai(0) = A 2 {0) = and A'^0) = A' 2 (0) = -1, 
Q implies that for all r G (0, T], 

Recall T < min{l, r }. It follows that for all r € (0, T], 

1 „ S 1,„ , O 1 q 1 



Mr) j 



Ai(r) 



< 9 ^ < * rf 9 ^ 2 + 

2 

By the definition of M, we have M < -M + \Sfi\, which implies 
(11) M<*\5n\. 



When 5[i = (i.e., /Ltj = /Lt 2 ), Inequality ( jlT] ) requires A/ = 0. Consequently, 
-^(t) = ^4i( r ) f° r an T G [Oj^ 1 ]- This proves the uniqueness of the solution on 
the interval [0, T]. Recall that with /ii G (/io — A*o + £)> ^i( T ) < for all 

tA' 

t G [T, To] by our choice of To and 5. Then A'" = with the initial conditions 

A(T) = Ai(T) < 0, A'(T) = A[(T), A"(T) = A'{(T) constitutes a regular initial 
value problem, to which A\ (r) is the unique solution throughout the interval [T, to] 
according to the standard existence and uniqueness theory. Altogether, we have 
shown the uniqueness of the solution on the full interval [0, To]. 

Since the choice of T is independent of n\ and /j, 2 , after combining (10 1 and 



(111 and integrating we can obtain a uniform constant Ct > such that for any 
Hi,H2 G (Mo -8,ho + S), 

\A 2 (t) - A x (t)| + \A' 2 (r) Ai(r)| + \A' 2 '{r) - A'{{r)\ < C T \ m M i| 

for all t G [0, T). Since Ai )2 (r) < q{r) < for all t G [T, to], by applying Gronwall's 
inequality on the interval [T, To] we can extend the above inequality to the full 
interval [0, To] after replacing Ct with some larger C > 0. □ 

Let Aq(t) be the solution to the IVP with A%(0) = fi G R. Clearly, we can 
extend Aq(t) to any f > To as long as Aq(t) < for all r G (0, f]. Furthermore, by 
tA' 

considering A'" — — — with the initial conditions specified at t and set close to 

OA . 

(Ao(to), A' (tq), Aq(to)), we can extend Lemma^to the interval [0, f]. In particular, 

by applying Gronwall's inequality on the interval [to,t], we obtain the following 

corollary. 

Corollary 3. For any f > such that Aq(t) < for all r G (0, f], there exist 
C > and a sufficiently small S > such that for any [i G (fxo — 5, //o + S) the 
solution A(t) with A"(0) = fi remains negative for all t G (0, f] and satisfies 

(12) \A(t) - A (r)\ + \A'(r) - A' (r)\ + \A"(r) - A '(t)| < C\n H 

for all t G [0, f] . 
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3.2. General behavior of solutions. Now that we know Q is well-posed, we 
demonstrate some properties of solutions when [i > 0. First, its solutions reach a 
first minimum in finite time. Let 

(13) n :=sup|r > : sup A'(s) < 1 . 

Note that T\ depends upon /x, though we supress this for notational simplicity. 

Lemma 4. For any /i > we have < t\ < oo. In addition, -A'(ti) = and 
A"(ti) > 0. Moreover there exist /zo,ao,ai > such that y, > fio implies (a) 

— < n < -, (b) -— < A(n) < -— and (c) /i < A"(n) < 2fi. 
2/j, /i A* M 

Proof. Clearly, -<4/(ti) = if tj is finite. For r € li := (0,Ti) (possibly n = oo), A 
and A' are negative and consequently A'" > 0. Thus A" is increasing. This implies 
that A"(t) > fj, > for r € Ji and also j4"(ti) > if ti is finite. Thus A is concave 
up on 1%. This together with the fact that A(0) — implies via the mean value 
theorem that \A(t)/t\ > \A'(t)\ for all t € h and so tA'/A < 1. Thus we have 

(14) < A" 1 < 1/6 

for t € I\. Integrating this inequality twice gives [it — 1 < A' < t 2 /12 + /.it — 1 
for all r £ Jj. If /i > 0, the lower bound on A' becomes positive in finite time, 
implying that t\ < oo. Examination of the zeroes of the upper and lower bounds 



gives the estimates for ri stated in the lemma. Further integration of ( 14 1 gives 
the estimates on A and A" . If fi = 0, the lower bound on A' is not sufficient to 
conclude that t\ < oo. However the fact that ^'"(O) = 1/6 can be used to improve 
the lower bound in this case. The details are standard and omitted. 

□ 

After reaching their first local minimum at r = Ti, solutions A either reach zero 
in finite time or alternately reach a local maximum (for which A < 0) in finite time, 
as we now demonstrate. Let 

(15) t 2 := sup \t > t 1 : inf A'(s) > and sup A(s) < 

^ T X <S<T T!<S<T 

As with n, T2 implicitly depends upon /i. If this is finite, then we have either 
A(t 2 ) = and A'{t 2 ) > or A(t 2 ) < and A'(t 2 ) = 0. 

Lemma 5. For all /i > we /icwe ti < t 2 < oo. 

Proof. Clearly r 2 > n since A"(ri) > 0. Suppose that r 2 = oo. Since A is bounded 
above and nondecreasing for r > t\, we have A(t) convergent as t — > oo. Since 
we have A < and A' > for all r > r 1; (7| implies A'" < for all r > r x 
as well. Thus A" is nonincreasing. Since it is initially positive and A converges, 
clearly A" must eventually become negative. Since A" is nonincreasing, once it is 
negative it remains negative for all subsequent times. Thus there is an e > and 
a value t% > T\ such that A"(t) < —e for all r > T3. Integrating this, however, 
indicates that A' must eventually become negative, a contradiction, and the lemma 
is shown. □ 
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3.3. Energy estimates and behavior of A for ji = or ji large. Our next goal 
is to show that for [i large, solutions go to zero in finite time but when fj, = the 
solution has a local negative maximum at r 2 . This requires the following energy 
estimates, similar to those for 

Lemma 6. Fix /i and take T\ and t 2 as above. Let 



Then for all t\ < r < t 2 



\[A'{t)Y KV^Air)). 



Moreover, take r :— 6A n '(ti) and set 



T 



V2(a) := — a In 



6 \A( n ) 

If f > r 2 then for all t± < r < r 2 

V 2 {A{t))<\[A\t)]\ 

Proof. Let J 2 := (t~i,t 2 ). By definition of r 2 , we have ^4(r) < and A'(t) > for 
all t € J 2 . If f > r 2 then we have for all t E I 2 , 

T*M < A »> (r) < ^(r) 
6^(r) " [Tj - 6A(t) ' 

where the equalities hold only if A 1 (r) = and the second inequality is true inde- 
pendent of f . Integrating this from t\ to r gives 

f - ^)<Aii(r)<A>i(ri) + T -LlJm) 



+ 6 H J < < + 6 ™ u(ri) ; 

for all t 6 7 2 . Note that we obtain strict inequalities here because with A"(t\) > 
there is a small a > such that A'(t) > for all r 6 (ti,ti + a). If we multiply 
by A'(r) and integrate once more from t\ to any r G (ti,t 2 ] then the conclusions 
of the lemma follow. In particular, f is irrelevant to the upper bound V±. □ 

Lemma 7. There exists M > sitc/i f/iaf fi > M implies A(t 2 ) = and v4'(t 2 ) > 0. 

Proof. First we show that if f > r 2 then we are done. We know by the definition 
of t 2 that either A'(r 2 ) = or A(t 2 ) — 0. Inspection of V 2 (a) demonstrates that 
V 2 (a) > for a g (A(ri), 0) and is zero at the ends of that interval (see Figure [3]). 
Thus, since A is increasing on (ti,t 2 ) we have A(t 2 ) € (^4(ri),0] and therefore 
V 2 (A(t 2 )) > 0. We also know for r g (ti,T2) that A' > 0. So in particular Lemma 
M implies A'(t 2 ) > ^2V 2 (A{t 2 )) > 0. So, since ^4'(t 2 ) ^ it must be the case that 
A.(t 2 ) = 0. 

On the other hand perhaps r 2 > f = 6A"(ti). We will prove by contradiction 
that this is impossible when /j, is large. The estimate in the proof of Lemma [6] gives: 

A'{t) 



1 < 
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Figure 3. Sketch of A'{t) vs A(r) when \i is large. 



which is valid for t± < t < r if we assume r < ti . (Note that for fi big enough, the 
estimates on t\ in Lemma [4] imply that f > n.) Integration from t% to f gives 

r A'(s)ds _ r A{f) da 
Tl - J T1 y/2V 2 (A(s)) ~ J A^) y/ZVoJa)' 
Substitution in from the definition of Vi and f gives 

1 r A & da 

v /2A"(r 1 ) A(rx) v/aln(a/A(ri)) : ' 
A change of variables and the fact that 1 > A(t)/A(t\) > implies: 



6A"(r 1 )-r 1 < 



6A"( n )-n < 
Letting X 



l^(n)| 



db 



< 



\A(n)\ 



db 



2A"( n ) Jaw/aw V\bMW\ ~ V 2A "( T i) Jo \/\bWfij\' 
db 

< oo we have 



o V\b~Hbj\ 



6A"(n) - n < 



l^(n)| 



_ 2A»(n)' 

Now we employ the estimates in Lemma [4] to find for all /i > /io: 

if /oo 
^ V 2 ' 
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A A' 



Figure 4. Sketch of A'(t) vs Afa when p = 0. 



For p sufhciently large this is impossible and so the lemma is shown. 
Lemma 8. If p = 0, then Afa) < and A'(t 2 ) = 0. 



□ 



Proof. Proving the lemma amounts to showing that Afa) ^ 0. Notice that Vi(0) = 
Afa) fa/6- A" fa)). Let pfa = t/6-A"(t). p(0) = 0, since p = 0. Also 
p'fa = 1/6 - A"'(t). In the proof of Lemma [4] we saw that A'" < 1/6 for r e /i 
(see fli| )) and so p' > 0. Consequently t x /& - A" fa) > 0. Thus, since Afa) < 0, 
we have Vi(0) < 0. Moreover, V x {Afa)) = and V{(Afa)) = A" fa) > 0. Thus 
there is a point k G (A(n), 0) at which Vi(k) = 0. See Figure|4j 

Since for r € [ri,r 2 ] we have > A(r) > A(r x ) and < A'(r) < ^2V 1 (A(t)) 
we must therefore have > k > Afa. Since t 2 is hnitc for p = and ^4 is bounded 
from zero for all t\ < r < t 2 , the definition (15) of t 2 requires ^4'(t 2 ) = 0. 



□ 

3.4. Openness of qualitative behavior. The next pair of results show that 
slightly varying p does not change the qualitative behavior of A at r 2 . That is 
to say, if po is such that Afa) — then there is an open neighborhood of po 
where the corresponding solutions enjoy this same property. Likewise, if po such 
that A' fa) = 0, then there is an open neighborhood of po where the corresponding 
solutions do the same. Let Afa) be the solution to the IVP ^ with A"(0) — p > 0, 
and define f 2 for Afa as in (15). 
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Lemma 9. Suppose A(t 2 ) = and A (r 2 ) = 2h > 0. Then for any e > there 
exists a 8 > such that for any fx g (fi — S, /t + 5), i/ie solution A(t) to the IVP Q 
im'tt A"(0) = /i reaches zero at t 2 < t% + e- 

Proof. Since A(t 2 ) = 0, the right hand side of Q is undefined at this point and so 
we cannot directly use our results about continuous dependence on /j, from Corollary 
[3]to prove this result. Instead we will have to rely instead on the nearly hamiltonian 
structure of Q. 

Integrating the ODE of Q with the initial conditions A(Q) = 0, A'(0) = -1, 
and A"(0) = /1, wc obtain that 



(16) 



A(r)A"(r) - \ (A'(r)f + \ = \tA{t) \ 



A(s) ds 



as long as A < on the interval (0, r). 

Since A(t2) = and A'(t2) = 2h > 0, we can choose a small a > such that 
\A(t)\ > \h(r - f 2 )| and A'(r) < 3h for all r g [f 2 - a,? 2 ]. It follows that 



A"(r) - A"(r 2 - a) 



< 



sA'(s) 



6A(s) 



ds < 



? 2 3h 



6h(s - f 2 ) 



A 1 - 1 Q 

as = -r 2 In ^ . 

2 t 2 - r 



Therefore, A(t)A"(t) —> as r — > r 2 . Applying (16) to A and taking the limits of 
both sides as r — > t 2 ~ , we obtain 



-2V 



A(s) as. 



Suppose there exist an e > and a sequence {/i n } — > A such that for any n, the 
solution A n {r) to the IVP ^ with A^(0) = fi n is below zero for all r g (0, r 2 + e). 
Since A(t) < for all t g (0,f 2 ), we can apply (12 1 of Corollary [i] to obtain a 
subsequence {Hn k } and some Q > such that for each k > Q, 

\A nk (r) - l(r)| + KM - A'(r)\ + \A"(r) - A"(r)\ < { 



for all r g [0,r 2 - i]. Clearly, A^(f 2 - i) -> A'(f 2 ) 
generality, we assume that A' (r 2 ) > /i for all k > Q. Let 



2ft. Without loss of 



sup { t g (f 2 - -,? 2 + e) 



inf 

1 

f 2 - - < s < r 



AL(s)>/i 



T 2 . 



Note that Wfc — >■ as k — >■ 00. Otherwise, for some very large k, A nk (r) would reach 

zero at some r £ (r 2 — — , t 2 + Wfc) since (r 2 — — ) converges to ^4(t 2 ) = from 

k ^ 
below. Without loss of generality, we assume that t 2 + cjf. < r 2 + e for all fc > Q. 

Then the definition of uik implies A' + oJk) = h an d A^ fc (T 2 + u)k) < for all 

k > Q. Applying ( 16 ) to A nk at r = r 2 + Wfc yields 



A„ fc (r 2 + w fe (f 2 + u k ) - ^ 2 + l - 



- (f 2 + w fe ) A„ fc (f 2 + UJ k ) - g 



^4n fc (s)ds. 
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Since > A nk (f 2 + u k ) > A nk (f 2 



A(t 2 ) — 0, the right-hand side converges 



to 



G 



A(s) ds 



-2h 2 



However, for all k > Q the left-hand side is greater 



1 2 1 

than or equal to — -h + — . This is a contradiction. 
In the next lemma, we still assume A"(0) = p > 0. 



□ 



Lemma 10. Suppose A(j 2 ) < and -A'(r 2 ) = 0. Then there exists a 5 > such 
that for any /i£ (p — 5, p + S) , £/ie solution A(t) to the IVP Q with A" {§) = fi 
satisfies A(t 2 ) < and A'{t 2 ) =0. 



Proof. Define n for A{t) as in ([13]). Then A"(?i) > since A"(0) = p > 0. It 
follows that A'(fi) = A'(t 2 ) = and A'(r) > for all r € (ri,f 2 ). Thus there 



is a t c € (ti,T2) such that A"(j c ) — 0. Furthermore, since A" 



6A 



< for 



all t e [t c ,t 2 ), we have ^4"(t 2 ) < 0. This allows us to choose a small a > such 
that A(r) < for all r e (0,f 2 + a], A'(r 2 + a) < 0, and A"(r 2 + a) < 0. Take 
a small j3 > 0. Then Corollary [3] guarantees the existence of a (5 > such that 
for any /j, £ (/t — <5, p, + S), the solution A(t) with A"(0) = /i is below zero for all 
r e (0,r 2 + a] and satisfies A'(r 2 — /3) > and A'(t 2 + a) < 0. This implies that 
A(t 2 ) < and A'{t 2 ) = at t 2 G (t 2 - /3, r 2 + a). □ 

3.5. Final steps. We are now in a position to complete the proof of Proposition [T] 

Proof. Define the set VL as follows: 

fi := {/2 G R : A(r 2 ) = and A'{t 2 ) > for any fi > p,} . 

Lemma [7] guarantees that f2 is nonempty, and Lemma [8] shows that zero is a lower 
bound of f2. Take /Lt* := inf > 0, with the inequality guaranteed by Lemma [8] 
and Lemm a [To] Let A*(t) be the solution to the IVP with A»(0) = /Lt». Then 
by Lemma [9 and the definition of /x*, it is only possible that either j4*(t 2 ) = and 
^4 (j 2 ) = or A*(r 2 ) < and A^(r 2 ) = 0. Furthermore, we rule out the second 
case by the definition of and Lemma 10 This proves (i)-(iv) of Proposition [T] 
Since A*(t) < and A' t (r) > for all r € (ri, r 2 ), we have 

t.K Tl A> 



6 A, ~'~* - 6 A, 
for all r € (ti,t 2 ). Integrating the above inequality gives 



n 

6 



In 



A,(r) 



A*(ti) 



+ <(r 1 )<<(r)<|ln 



+ <(ri) 



for all r G (ti, t 2 ). Thus A"(r) — > —00 as t — > t 2 . Then we can take a sufficiently 
small a > such that A"(t) < — 1 for all r G (r 2 — a,r 2 ). By the mean value 
theorem we have 

K(t 2 ) - AUt) = A"(s)(t 2 -t)< -(t 2 - r) 
for any r £ (r 2 — a, r 2 ) and some s € (r, t 2 ). Since A^(r 2 ) = 0, this implies that 
< (t 2 - t) < A'^(t) and consequently -(r 2 - r) 2 < A*(r 2 ) - A*(t) = -A*(r) for 
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all r e (r 2 - a,r 2 ). Thus ln|A»(V)| = 0{\n(\r - r 2 |)) as t 
of Proposition [T] 



r 2 . This proves (v) 



□ 



4. Numerical Study of K(2, 2) 

4.1. Regularization and scaling. In this section we numerically assess the ill- 
posedness of pj). But it is inappropriate to directly simulate an equation that not 
only lacks a local well-posedness theory but for which we suspect ill-posedness. Thus 
we regularize the equation. Simulating this regularized problem, we find evidence 
of the ill-posedness as we let the regularization parameter vanish. 
We study the following regularization of ([2| 

(17) (7 + Sd 4 x )d t u = d x (u 2 ) + &l (u 2 ) . 

Here I is the identity and 6 > is the small regularization parameter. We choose 
this particular regularization since its implementation is natural when simulating 
solutions of Q via a pseudospectral method. 

Inverting the operator on the left-hand side to be able to write it as an evolution 
equation, we see that d x /(I+6d*) and 9^/{I+6d£) are bounded operators. Indeed, 



< CVT 



and 



I + Sdi 



< CVT 



The boundedness of these operators make it trivial to prove: 

Theorem 11 (Local Well-Posedness of a Regularized Problem). Fix S > 0. Then 



(17 1 is locally well-posed in H s for any s > \ 

• For all uq € H s , there exists a T > and a unique function u{t) € 



C(0,T; H s ) solving the integral form of (17), 

• The solution will depend continuously upon the data, 

• There exists a maximal time of existence, T cx i st such that if it is finite, 

lim 



Ht)\\ H s 



oo 



The proof of this, which we omit, is an elementary application of the fixed 
point technique. The operator is bounded in L 2 -based Sobolev spaces. For s > \, 
H s is an algebra which makes the nonlinearity easy to treat. An impediment to 
extending Theorem [TT] to one which is global in time, or even one which holds on 
time intervals of a length uniform in 6, is the lack of an obvious coercive conserved 



quantity associated with (17). The only obvious conserved quantity for (17) is 



u(x,t)dx = / u(x,0)dx. 



This is formally an invariant for solutions of (|2|) as well 



By simulating the regularized equation ( |17[ ) we seek evidence of ill-posedness of 
([2]). In particular, we shall find a sequence of vanishing initial conditions for (17), 
such that as 8 — > 0, the corresponding solutions at t = > have H 2 norm of 



(at least) unit size. To that end, we will scale the initial data for (17) using the 
same scaling given by ([5| and link the vanishing of the smoothing parameter S to 
the scaling parameter A. Specifically we will study (17) with 



(18) 



S = 5\ := .IX and 



l (x,0) = h(x) :=A/(^) 
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where 

f{x) = exp(-4a; 2 ) 

and A — > 0. Note that the calculations leading to ^ show that H/xllip < CA 7 / 6 
and so this choice of initial data vanishes in the limit. We denote the solution of 



(17 1 with the choices in (18) by u\{x, t). Note that we have taken our initial data 
to be everywhere positive, unlike the self-similar solutions found for p}. We do 
this to demonstrate that the problematic effects of degeneracy manifest themselves 
even for solutions which do not cross the a;-axis. 

The choices for 8\ and fx (x) are made for the following reason. Instead of scaling 
just the initial data for (17) suppose that we rescale the whole equation by 



(19) 



Then (17) becomes: 



ux(x,t) = Xv x (tt73>*) 



where y — x/X 1 ^ 3 . Plugging in from (18) we get 



vx(v,0) = f(y) 



and 
(20) 



(I- 



AX^d^vx 



X 2 ' 3 d y v 



Of course (20) is not an exact scaling. Notice that as A 



the term d y v 2 term 
will vanish. The equation will asymptotically be dominated by the third derivative 
term which is consistent with our expectation that this term is the source of the ill— 
posedness. Moreover notice that the coefficient of regularization even in these scaled 
coordinates vanishes. That is to say, we have chosen the regularization parameter 
so that it vanishes "more rapidly" than scaling effects. 

4.2. Computational results. We integrate our problem pseudospectrally with a 
Crank-Nicholson time stepping scheme. Some additional details are presented in 



Section|473] Simulating ([17} with Q to T = .1 on the domain [-2rr, 2rr) with 8192 
resolved Fourier modes, the data evolves as in Figures [3] and [6j At large values of 
A, there is some development of oscillations. At very small values of A, it develops 
a highly oscillatory structure. Note that these oscillations appear for x > 0, where 
/' < 0. That is to say in the exactly the place where the term 6u x u xx in |2]) acts, 
heuristically, as a backwards heat operator. 

We present the evolution of the L 2 , H 1 , and H 2 norms in Figure[7|of the solutions 
as functions of time. Though there is little growth in time of the L 2 norm, we see 
orders of magnitude jumps in the H 1 and H 2 norms at t — .1 as we send A — » 0. 
We contend that this is evidence of ill-posedness. We have a sequence of initial 
conditions fx, which have vanishing H 2 norm. At a fixed T > 0, the H 2 norms 
are growing as we let A — > 0. Thus, there is a loss of continuous dependence of the 
solution map upon the initial data about the u{x, t) = solution. 

This growth in the H 2 norm corresponds to a spreading in Fourier space. Indeed, 
Figure [8] shows precisely this. At fixed time, the Fourier support grows as S — >• 0. 
These figures also indicate that except for A = .04, all simulations are extremely 
well resolved. Even the A = .04 simulation is well resolved through t = .06, and 
only slightly under resolved at the final time, t = .1. 




4.3. Details of numerical methods. As noted, the simulations presented in Sec- 
tion [4] were generated using a fully de-aliased pseudospectral discretization with 
Crank-Nicholson time stepping. This was performed in Matlab and f solve was 
used to solve the nonlinear system at each time step with a tolerance of 10 -8 . 8192 
modes were resolved in these simulations, with a time step of At = .001. 
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0.1 - 



0.1 - 



: 0.05 




0.05 






x 

(a) 
( = 0.04 



0.5 -0.5 



(b) 



0.1 - 



0.1 - 



3 0.05 






x 

(c) 

t = 0.08 



S 0.05 




0.5 -0.5 



(d) 



0.1 - 



0.05 




0.1 - 



S 0.05 




(c) 



(f) 



Figure 6. Evolution of the regularized K(2, 2) equation with data 



given by (18) with A = .05. 



The invariant / u(x,t)dx was conserved to at least eight significant figures in 
all the computed cases, as shown in Table [T] 

4.4. Convergence of Norms. Since we are concerned with the value of H 2 at 
the end of our simulation, it is also important to confirm that this is converging to 
a fixed value as we refine the spatial and temporal resolution of our simulations. 
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Figure 7. Evolution of several norms for the regularized K(2, 2) 
equation with data given by (18). 

Table 1. Numerical approximation of J u(x,t)dx at the out- 
putted times for different values of A 



0.00 0.886226925453 

0.02 0.886226925477 

0.04 0.886226925489 

0.06 0.886226925478 

0.08 0.886226925497 

0.10 0.886226925509 



0.4 

0.261191032665 
0.261191032667 
0.261191032674 
0.261191032676 
0.261191032673 
0.261191032673 



A 
0.2 

0.103653730004 
0.103653730004 
0.103653730004 
0.103653730004 
0.103653730004 
0.103653730004 



0.1 

0.0411350100123 
0.0411350100123 
0.0411350100123 
0.0411350100123 
0.0411350100123 
0.0411350100123 



0.05 

0.0163244395416 
0.0163244395416 
0.0163244395416 
0.0163244395416 
0.0163244395416 
0.0163244395416 



For the initial condition with A = 0.2, Table [2] gives the values of the H 2 norm for 
various choices of the number of resolved Fourier modes, N, and the time step, At. 
With only 512 resolved modes, we appear to be fully resolved in space, and there is 
little gain in accuracy when N is increased. Of course, for other initial conditions, 
more than 512 modes may be needed to resolve the simulation. 

For the three values of At, we appear to have achieved at least one significant 
digit of accuracy and as we reduce the time step, the variations of H 2 become 
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= 0.1 


• A 
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(c) 



(f) 



Figure 8. Evolution of the regularized K(2, 2) equation with data 
given by (18) with A = .05. 



smaller. Since we are using a Crank-Nicholson scheme, we expect 0(At 2 ) conver- 
gence. If we perform Richardson extrapolation, the estimated value of H 2 based 
on the At = .004 and At = .002 simulations is 697.87. Comparing the At = .002 
and At = .001 simulations, the estimated value is 697.97. Separate computations, 
performed with just 512 grid points for expediency, show that with At = 0.0005 
the norm takes the value 697.96, and for At = 0.00025, the value is 697.99. 
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Table 2. Values of H 2 at t — .1 corresponding to the A = 0.2 
simulation with different numbers of resolved Fourier modes, N, 
and time steps, At. 







At 






0.004 


0.002 


0.001 


512 


675.27623578 


692.22147559 


696.5301160 


1024 


675.27623585 


692.22147576 


696.5301162 


N 2048 


675.27623587 


692.22147576 


696.5301162 


4096 


675.27623586 


692.22147576 


696.5301162 


8192 


675.27623586 


692.22147576 


696.5301162 



Table 3. Values of H 2 at t = .1 corresponding to the A = 0.1 
simulation with different numbers of resolved Fourier modes, N, 
and time steps, At. 







At 






0.004 


0.002 


0.001 


2048 


12452.791732 


14474.328299 


15077.20691 


N 4096 


12452.791732 


14474.328298 


15077.20691 


8192 


12452.791732 


14474.328299 


15077.20691 



An examination of the A = 0.1 case is similar. For the same three values of 
At, the H 2 norms, given in Table [3J are in agreement on the order of magnitude. 
Richardson extrapolation using the data at At = .004 and At = .002 predicts a 
value of 15148, while the predicted value based on the At = .002 and At = .001 
value is 15278. Independent simulations with 2048 grid points yield a value of 15253 
when At = 0.0005 and a value of 15302 when At = 0.00025. Thus, we believe that 
our choice of discretization parameters, N = 8192 and At = 0.001, for the results 
presented in Section 4.2 yield meaningful measurements of H 2 that are accurate to 
at least one signfigant figure. This is sufficient for our study of ill-posedness. 
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